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SYSTEM FOR PHASE ERROR ESTIMATION FOR A QUADRATURE MODULATOR 

FIELD OF THE INVENTION 
The present invention relates generally to error estimation systems and methods and more 
specifically to a method and system for estimating phase error inherent in waveforms sent from 
analog transmitters. 

BACKGROUND OF THE INVENTION 
Transmitters known to the art are driven by digital circuits at very high sampling rates. With 
these digital circuits, the inphase and quadrature data may be created digitally and in perfect 
quadrature. However, the inphase and quadrature data may be distorted when utilized with analog 
circuits during transmission due to phase, gain, and direct current error inherent in such analog 
circuits. 

A compensation of the phase, gain, and direct current error is necessary for clear reception of 
the transmitted signal. Another problem with estimating phase errors lies in the variation of the 
error. Phase error may change and may be altered due to variations over time, temperature, and 
wider bandwidths. 

Consequently, it would be advantageous if a system and method existed that could estimate 
the phase error in the waveform sent from an analog transmitter. It would also be advantageous if the 
system and method could continuously monitor the phase error due to variations in phase error. 

SUMMARY OF THE INVENTION 

Accordingly, the present invention is directed to a novel system and method for continuously 
estimating the phase error inherent in a waveform sent from an analog transmitter. Alternative 
approaches may be utilized to estimate the phase error. A first approach involves computing the 
minimum and maximum envelope values and computing a ratio of the minimum and maximum 
envelope values. A relationship between the ratio of minimum to maximum envelope value data and 
phase error may be determined. Other alternative approaches involve utilizing a half angle technique 
and a maximum abscissa technique for estimating the phase error. 

It is to be understood that both the foregoing general description and the following detailed 
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description are exemplary and explanatory only and are not restrictive of the invention claimed. The 
accompanying drawings, which are incorporated in and constitute a part of the specification, 
illustrate an embodiment of the invention and together with the general description, serve to explain 
the principles of the invention. 

BRIEF DESCRIPTION OF THE DRAWINGS 

The numerous objects and advantages of the present invention may be better understood by 
those skilled in the art by reference to the accompanying figures in which: 

FIG. 1 depicts an embodiment of a flow diagram for estimating phase error in a waveform by 
the minimum to maximum envelope ratio technique of the present invention; 

FIG. 2 depicts an exemplary table of values for comparison; 

FIG. 3 depicts an exemplary graph for determining the phase error based upon a given 
minimum to maximum envelope ratio; 

FIG. 4 depicts an embodiment of the distortion after correction utilizing a straight line 
method and a quadratic equation; 

FIG, 5 depicts an embodiment of a flow diagram for estimating phase error in a waveform by 
the maximum abscissa technique of the present invention; 

FIG. 6 depicts an exemplary graph for determining the phase error via the maximum abscissa 
technique of the present invention; 

FIG. 7 depicts an embodiment of a flow diagram for estimating phase error in a waveform by 
the half angle technique of the present invention; and 

FIG. 8 depicts an exemplary graph for determining th$ phase error via the half angle 
technique of the present invention. 

DETAILED DESCRIPTION OF THE INVENTION 

Reference will now be made in detail to a presently preferred embodiment of the invention, 
an example of which is illustrated in the accompanying drawings. 

Referring to FIG. 1, an embodiment of a process 100 of the present invention is shown. The 
process 100 may begin by continuously determining the minimum and maximum envelope values 
for the waveform 1 10. The envelope of a waveform refers to a point-by-point magnitude of the 
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waveform. A ratio of the minimum envelope value and the maximum envelope value may be 
computed 120. Utilizing the ratio of minimum envelope value to the maximum envelope value, a 
comparison between the ratio and the phase error 130 may be made. The relationship between the 
ratio and the phase error can be computed and data may be fit to a straight line or, preferably, a 
quadratic equation. Fitting the data to a quadratic equation may yield a more accurate estimate of 
the phase error. An advantage of the process 100 is ability to not require quadrature components. 
This is advantageous as less computations may be required in the estimation of the phase error. 

A simulation of the method of estimating phase error by utilizing the ratio of minimum 
envelope value to the maximum envelope value may be accomplished utilizing the following 
exemplary factors. 

A number N data points may be created for study: 

N:=1000 

n :=0..N- 1 

An angle increment may be created: 

A '= 360 . n 
N 180 

9 :=n-A 
n 

The angle may be converted from radians to degrees by 

The semi major and minor ellipse axes and direct current offsets may be: 

a:=l 

b:=l 

h:=0 

k:=0 

The x-axis variables may be: 
x mn :=-a+h 

x mx :=a+h 



5 x :=2-a-N 
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1 



V =x mn+ n ' 8 x 

A number of phase errors up to 20 degrees may be created: 
<tmax:= 20 

nl := 0..(|tnax 
♦hi 1 

A matrix of waveform coordinates, u and v, may be created with a column of N points for each phase 
error (the phase error may be inserted in the v channel): 

u t := a- cos (9 n ) 
n,nl v "/ 

V n,nl : " b ' sin ( On + ^0 

The envelope of the coordinates may be computed by: 

env , := Hu ^ 2 + (v t ) 2 
n,nl ^\ n,ni; \ n,nlj 

The maximum envelope value for each phase error may be computed by: 

max(env^ n ^) 



emx . 
nl 



The minimum envelope value may be computed by: 
:= min(env^ n ^) 



emn , 
nl 



The ratio of minimum to maximum envelope values may be determined by: 



emn . 
nl 

nl emx 



nl 

The phase at the maximum envelope value may be calculated by differentiating the expression for 
the envelope squared with respect to 0, setting it to zero, and solving for 0: 
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( 2 i x ^ 
b .sin^.^!) 

atan — 

^a -b -008(2-^1) j 
^nl : ~ ~ 

The x and y coordinates of the maximum envelope at 0 may be determined by: 
X nl := a.cos(0 nl ) 

Y ni :=b,sin (©nl + <l>nl) 

The x and y coordinates of the minimum envelope at 0 may be determined by: 
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Xmn nl :=a*cos| Q n \ + — 



Ymn t := b-sin 
nl 



®nl + " + <t>nl 



The maximum envelope (envMx) and minimum envelope (envMn) may be computed by: 
envMx nr=J( X nl) 2 + ( Y nl) 2 

envMn nl :=J(Xmn nl ) 2 + (YmnJ 2 

The maximum and minimum envelope values so computed (envMx and envMn) may be 
compared with the minimum and maximum envelope values computed previously (emx and emn). 
Referring to FIG, 2, an exemplary table showing a comparison of the separate values for the 
minimum and maximum envelope values is shown. The values computed in the exemplary table 
may be calculated as follows: 
R nl)0 :=^l-C 

R - , := emx , 
nl,l nl 

R , „ := emn , 
nl,2 nl 



R , * := er nY 



nl 



-5- 



PATENT 
00CR020/KE 



R 



nl,4 



R 



nl-,5 



:~ envMx 



nl 



nl,6 



:= envMn 



nl 



The values computed for envMx and envMn are the same as the values for emx and emn, showing 
that the equation for 0 is correct. 

In order to estimate a phase error from a minimum to maximum envelope ratio, a 
relationship between the ratio and phase error may be computed. Utilizing a curve fitting technique, 
data may be fit into a straight line and more preferably, a quadratic equation. 

Beginning with a straight line curve-fitting technique utilizing the R n i 3 values from FIG. 2, 
a rise per degree of run, slope of ratio vs. phase error in degrees may be shown by: 



mrd= -0.014989570215612 

The straight line fit of ratio may be shown by: 

intcept := 0.9945 

er est := mrdcLi-C + intcept 

nl 

The above equation may be used by solving it for § and substituting the computed ratio value 
from the measured minimum and maximum values of the envelope (emx and emn), The resulting § 
is then the phase error estimate, <|> es t. 

(()est = (er est - intcept) -mrd 1 

The phase error may be computed in terms of the computed minimum to maximum envelope 
ratio. The computed minimum to maximum values with the related ones in the table as shown in 



mrd:= 




<fmax 
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FIG. 2 along with the exemplary plot as shown in FIG. 3 may show the error in the estimate. For 
plotting purposes the computed envelope ratios may be used as measured values. 



<Mnl := 



ernx , - intcept 

nl 

mrd 



Rise per radian of run, slope of ratio vs. phase error in radians: 



er nXj 



mr:= ■ 



(frmax 



- 1 



(|max^ 

mr= -0.858839110069583 



er n x nl - intcept 

(jestn! := C 

mr 





il-o 














7 


.8 






; , f ;.^;., : 


P 


-0.367 


0.787 


1.922 


3.038 


4.135 


5.215 


6.277 


7.323 


8.353 


9.368 


10.367 


11.352 



T 

([est = 



The equation for (|) es t can be used as described above to estimate the phase error from a 

measured and computed ratio, or else the plot below can be used as follows. Draw a horizontal line 
at the computed ratio and find its intersection with the line used to fit the data. Then read the phase 
error value at that intersection from the abscissa. Referring to FIG. 3, if ratio = 0.9, the x coordinate 
at the intersection with the line for er es t is 6.277 where 6 is the exact answer. 

Assuming perfect gain and direct current offset correction so that only residual phase error 
remains and assuming straight-line curve-fitting the distortion caused by residual phase error can be 
computed using equation 8.40 from: 

Rockwell Corporation, Collins Division, Single Sideband Systems and Circuits 2 nd ed. 
(McGraw-Hill, Inc. 1995) 

distl nl + 1 " 2 * C0S (<MnK~ 1 " <kil)Xl + 1 + 2.cos((|est n K~ 1 - 4w)) 
stLineDistDb nl := if^distl nl = 0,-65, 10 log^distl^ 
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maj<stLineDistDb) = -49.875 

Hence, in one iteration, the phase estimate may be corrected to about -50 db using simple 
straight-line curve-fitting. This has been done above by changing the y intercept so as to allow 
bipolar error where the zero error point may be obtained somewhere between the endpoints of the 
curve instead of exactly at the endpoints. Distortion may be expected to vary from -50 db to -70 db 

with one iteration of this technique as the phase error varies from a magnitude of 0° up to 20°. 

Since residual gain and direct current offset errors may degrade the phase estimate, it is presently 

preferred to limit the estimates in practice to remain on the fitted curve. Again, a quadratic equation 

may provide a better fit and may provide better distortion characteristics. As shown next, the 

improvement in performance may be more than 20 db using a quadratic equation. 

Let the x-axis data be: 
xx:= <|>-£ 



Fit to the y data: 

( 



G(z,B) 



B A + B,-z+ B„-z 



"0 



vg2:= 



1 



B := genfit(xx, er nx , vg2 , G) 



( 0.999573323 
-0.017148083 



V1.09964456& 10 



.-4 



gl nl :=B 0 + B r xx 1 + B 2 .(xx nl ) 2 
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-B, - 




nl 




- 1 



1 + 1 + 2-C0S Xe St < - fei 



ma?(quadDistDb) = -71.495 

Referring to FIG. 4, an exemplary graph describing the distortion after correction versus the 
phase error is shown. As shown in FIG. 4, utilization of the quadratic equation of the present 
invention provides more accurate results. 

An alternate method of estimating the phase error in a signal from an analog transmitter in 
accordance with the present invention is a maximum abscissa method. Transmitter and receiver 
waveforms are elliptically distorted by phase and gain errors. The line drawn from the origin 
subtending angle (phase error) with the x-axis may intersect the distorted waveform at or near the 
extreme x value for phase errors less than or equal to 20°. The phase angle, <|>, at which the 
intersection may occur is derived by setting v = w (see equations below) where x = u at the point of 
intersection on the ellipse and solving for <[). 

Assuming direct current offset correction has already occurred this technique may use the 
point of tangency of a vertical line at x = max(I) to give the I and Q coordinates that approximate the 

phase error (atan[Q/I]) for errors of magnitude less than or equal to 20°. Sequential iterations 
converge on the actual phase error. This technique may have the advantage of being quick as the 
maximum x value may be found and the angle subtended by a line from the origin to this point on 
the waveform may be determined. The disadvantage is that it is only useful for small angles less 
than about 20°. 

Referring now to FIG. 5, an embodiment 500 of the maximum abscissa technique for 
estimating phase error is shown. The technique may begin by determining the maximum abscissa 
value for the waveform 510. The index of the maximum abscissa value may then be computed 520. 
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Once the index has been computed, the phase error may be estimated by finding the arctangent of an 
angle formed by a line from the origin through a point on the waveform at the index 530. 

A simulation describing the process may begin by selecting a phase error in degrees (if 1 
converts it to radians): 

N data points may be created for study: 
N := 1000 
n :=0.. N- 1 

An angle increment may be created by: 

N 180 
9 n :=n.A 

The angle increment may be converted from radians to degrees by: 
? := lSOn " l 

The semi major and minor ellipse axes and direct current offsets may be created by: 

a :=1 

b :=1 

h:=0 

k:=0 

An x-axis variable may be created by: 
x mn :=-a-hh 

x mx :=a-hh 
5 X :=2*a-N~ 1 
V =x mn+ n ' 6 x 

A reference line at a 45° angle may be created for plots below by: 
y45 :=x 

J n n 

A line from the origin at angle, <|>, and slope tan(<|)) may be created by: 

10 - 
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w n :=tan((j> )*x n 



The x (or I) component of distorted waveform may be computed by: 
u n :=a.cos(9 n )H-h 

The y (or Q) component of distorted waveform may be computed by: 



v n :=b.sin(e n -H)-hk 



©, the angle @ in degrees where v = w and the line w may intersects the ellipse waveform may be 
shown by: 



X may be on the line w if Y and Z are the same: 
Z :=tan((j) ) X 

which may show that w - v at 0 as desired. 
Y-Z = 0 

The maximum abscissa value may be found by: 

ux:=ma?(u) 

ux= 1 

The index of the maximum abscissa value may be found by: 
undxv n 1= ifj^u^ux, n, - 1 j 

undx := mso{ undxv ) 
undx = 0 

The phase estimate may be made based on the arctangent of the angle formed by a line from 
the origin through the point on the waveform at index undx. 




The X and Y coordinates at intersection may be given by: 
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♦ est '5 = 9.85107611658391 

While assuming perfect gain and direct current offset correction so that only residual phase 
error remains, the distortion caused by residual phase error can be computed using equation 8.40 
from: 

Rockwell Corporation, Collins Division, Single Sideband Systems and Circuits 2 nd ed. 
(McGraw-Hill, Inc. 1995) 



distDb := lO log(dist) 
distDb =-57.724 

Hence, in one iteration, the phase estimate for a phase error of 10° as in the above example 
may be estimated to approximately -58 db total distortion. Distortion may be expected to vary from 
-40 db to -1 1 7 db with one iteration of the maximum abscissa technique as the phase error varies 

from a magnitude of 20° down to 1°. Further iterations may improve the estimate, however, the 
error in the estimate could also be characterized by lookup table or the data may be fit to a quadratic 
equation that could be used to correct the estimate. Residual gain and direct current offset errors 
degrade the phase estimate which may be computed utilizing the maximum abscissa technique. 

Another technique for estimating phase error inherent in a waveform sent from an analog 
transmitter may be the half angle by intersection using envelope error technique of the present 
invention. Utilizing an envelope error may allow computing a phase error estimate, (|) es t, as twice 
the angle subtended by a line from the origin to the point of intersection of the received envelope 
with a circle of magnitude 1. The interpolation term may be added or subtracted depending upon the 
direction of rotation and the quadrant of the intersection. An advantage of the half angle technique 
lies in the accuracy of the estimated phase error. Another advantage of the half angle technique may 
be in the computation of an accurate estimated phase error even with direct current offset in the y or 
Q channel. 

Referring now to FIG. 7, an embodiment 700 of half angle by intersection technique for 
estimating phase error inherent in a waveform is shown. The technique may begin by calculating the 
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envelope for the waveform 710. The envelope error may be computed 720 and the indices of 
intersection may be found 730. Utilizing the indices of intersection and envelope error at sampling 
points on envelope on either side, an interpolation of an abscissa and an ordinate may be performed 
740. The angle subtended by a line from the origin to an intersection and the closest axis may be 
determined where the phase error is twice the angle formed 750. 

A simulation describing the half angle technique may begin by selecting N data points for 

study: 
N := 1000 
n:=0..N- 1 

An angle increment may be created by: 

A ■= 36Q 71 
N 180 

n 

The angle may be converted from radians to degrees by: 
C -180-tT 1 

The semi major and minor ellipse axes and direct current offsets may be shown by: 

a:=l 

b :=1 

h:=0 

k:=0 

An x-axis variable may be created by: 
x mn :=-a-hh 

x mx :=a+h 
V =x mn+ n5 x 

A phase error in degrees may be selected, in this simulation 10° is selected converts it to 
radians): 
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+ :=i<K 

The x (or I) component of distorted waveform may be computed by: 



u n :=a-cos 9 n th 



The y (or Q) component of distorted waveform may be computed by: 
v n :=b-sin(e n -H )-hk 

The envelope of u and v may be created by: 



4v 



The undistorted envelope would have a value of 1 .0. The envelope error may be computed by: 



e e :=en - 1 



The indices of intersection may be found by: 



index := 



*r0 

flag<-l 

ex<-0 

for v e 0.. N- 1 



if [ <v * i) -(%"\-, <o )] 



ex«-v 
i 



w-i-f 1 
flag<-0 



ex 



237 
487 

index = 

737 
987_ 

n2 :=0.. Iast( index) 

Using the indices of intersection and envelope error at sampling points on the envelope on 
either side of intersection to improve estimate via interpolation gives accurate estimates for angles 

from -90° to +90°. The quadrant and rotation may be solved by adding if envelope error goes from 
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positive to negative across the intersection, while subtracting if the envelope error goes from 
negative to positive across the intersection. The half angle technique may not produce an accurate 
phase error estimation with direct current offset in both channels in current form, but the output of 
the half angle technique could be ignored until most or all of direct current errors have been 
removed. If direct current offset exists in only one channel, then averaging all the <|> es t may provide 

accurate results as before without averaging and without direct current offset. 

The total change in envelope error between the two samples straddling the intersection of the 

distorted waveform with the reference may be computed by: 



A n2 := 



e e/mdex 



n2 



-emdex n2 -l 



The change in abscissa between the two samples may be computed by: 
The change in ordinate between the two samples may be computed by: 



*y B 2 :=v . 



index 



v/. 



n2 



index^-1 



The abscissa may be interpolated just before the intersection to provide an accurate coordinate at the 
intersection: 



for ke 0.. last( index) 



^"(inde^-l^- 



Andex^— 1 



■H if [ e e M <0] 



^-"(inde^-l) 



(inde^-1) 



•85^ otherwise 



x3 



The ordinate may be interpolated just before the intersection to provide an accurate coordinate at the 
intersection: 
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for ke 0.. last( index) 



y 3 k«- v («fa k -i)+- 



y 3 k < - V ^ndex k -l 



'efindex, -1] 



5 ^index^— 1 



•Sy 



^ if hN <0 ] 



■•8y. otherwise 



y3 

The angle of line from the origin to the point of intersection may be computed by: 



<b2 A i=atan 
e n2 



y2 



is 



n2 



x2 



is 



n2 



The phase estimate may be twice the angle subtended by a line from the origin to the intersection 
and the closest axis as shown by: 



♦ 2 est 



for ke 0.. last( index) 

♦ 3^2.(0.5., - *2 e J if [(♦2, t >°) (] y2 \ H 



♦ V-2--0.5-K-42, 



* 2 e<° V 



y2 



x2 



4> 3, — 4» 2 e -2 otherwise 



*3 



* 2 est^ = 



10.000037 
10.000039 
10.000037 
10.000039 



All of the estimates may be averaged to provide a correct result if there is direct current offset on 
only one channel: 



avg_*2 est :=C-(last(index)+ l)"'-/>J* 2 est 

\n2 / 



avg_<[)2 est = 10.0000382038299 

By assuming perfect gain and direct current offset correction so that only residual phase error 
remains the distortion caused by residual phase error may be computed using equation 8.40 from: 
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Rockwell Corporation, Collins Division, Single Sideband Systems and Circuits 2 nd ed. 
(McGraw-Hill, Inc.1995) 

dist := (l + 1 - 2-cos (avg_<|> 2 est 1 - + ) ) • (l + 1 + 2-cos (avg_<i> 2 est ^' 1 - + ) ) 

distDb :=101og(dist) 
distDb = - 129.542 

Hence, in one iteration, the phase estimate for a phase error of 10° as in the above example 
may be accurate to -130 db total distortion. Distortion may be expected to vary from -1 10 db to -135 
db with one iteration of this algorithm as the phase error varies from a magnitude of 20° down to 1°. 
Residual gain and direct current offset errors may degrade the phase estimate as computed by the 
half angle technique unless, as stated, there is direct current in only one channel (y or Q axis). 
However, the half angle technique may still be used if there is direct current in the x axis by shifting 
the origin by the amount of the offset in the x direction. 

It is believed that the present invention and many of its attendant advantages will be 
understood by the foregoing description, and it will be apparent that various changes may be made in 
the form, construction, and arrangement of the components thereof without departing from the scope 
and spirit of the invention or without sacrificing all of its material advantages. The form herein 
before described being merely an explanatory embodiment thereof, it is the intention of the following 
claims to encompass and include such changes. 



